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We develop a theory ol the tunneling spectroscopy for superconducting topological insulators 
(STIs), where the surface Andreev bound states (SABSs) appear as helical Majorana fermions. 
Based on the symmetry and topological nature ol parent topological insulators, we find that the 
SABSs in the STIs have a profound structural transition in the energy dispersions. The transition 
results in a variety ol Majorana fermions, by tuning the chemical potential and the effective mass 
ol the energy band. We clarity that Majorana fermions in the vicinity ol the transitions give rise to 
robust zero bias peaks in the tunneling conductance between normal metal/STI junctions. 
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Topological superconductors (TSCs) are a new state 
of matter [l]-Q that characterized by non-zero topolog- 
ical numbers of the bulk wave functions. They support 
topologically protected gapless surface Andreev bound 
states (SABSs), and the superconductivity infers that the 
gapless SABSs are their own antiparticles, thus Majo- 
rana fermions The realization of TSC and Majorana 
fermions in condensed matter physics is of particular in- 
terest because of their novelty as well as the possible 
application for quantum devices [H-[2lJ. 

The recent discovered superconductor Cu.j;Bi2Se3 
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25| is an intriguing candidate of the TSC because it is 
associated with another new state of matter, topologi- 
cal insulator: The parent material Bi2Se3 is originally 
a topological insulator with topologically protected gap- 
less Dirac fermions on its surface. With intercalating Cu, 
the superconductivity appears. On the theoretical side, 
it has been expected that Cu x Bi 2 Sc3 is a TSC by the 
Fermi surface criterion 26-28J, and possible SABSs spe- 
cific to this material have been studied Recently, 
a point contact spectroscopy experiment on this mate- 
rial has been done, and reported a pronounced zero-bias 
conductance peak (ZBCP) |30j]. With careful analysis 
excluding other mechanisms, it has been concluded that 
the ZBCP is intrinsic and signifies unconventional super- 
conductivity [3(|. Moreover, similar ZBCPs have been 
observed by other groups independently 33-351]. 

Motivated by this finding, we develop in this letter 
a general theory of Majorana fermions in superconduct- 
ing topological insulators (STIs) and their relation to the 
tunneling conductance. Up to this time, the relation be- 
tween SABSs and the tunneling conductance has been 
well understood in quasi two-dimensional superconduc- 
tors [l|: (1) If the SABS has a flat band dispersion as a 
function of the momentum parallel to the surface, k y , the 
resulting line shape of conductance always has a sharp 
ZBCP as realized in high-T c cuprate [J Hfj|. (2) If the 
SABS has a linear dispersion as a function of k v , the re- 
sulting line shape of conductance has a broad peak as ob- 



served in Sr2Ru04 37-33]. On the other hand, in three- 
dimensional superconductors, little is known about the 
relation between SABSs and the tunneling conductance. 
The only exception is a study on the superconducting 
analog of superfluid 3 He B phase Like Cu2;Bi2Se3, it 
is a three dimensional TSC supporting helical Majorana 
fermions on its surface |4l| - t46| . However, the resulting 
tunneling conductance always shows a double-peak struc- 
ture and never shows a ZBCP [13]. Therefore, in order 
to pursue the origin of the observed ZBCPs in STIs, one 
needs to develop a theory of the tunneling conductance 
for STIs. 

In this letter, we study the tunneling spectroscopy and 
underlying SABSs in STIs. Based on symmetry and 
topological nature of parent topological insulators, it is 
shown that SABSs in STIs generally have a profound 
structural transition of the energy dispersion (FigEJ. 
The transition results in a variety of helical Majorana 
fermions in SABSs, which we call cone, caldera, ridge and 
valley. We clarify that the transition naturally explains 
robustness of the ZBCP in STIs. From explicit calcu- 
lation, it is found that the tunneling conductances be- 
tween normal metal/STI junctions ubiquitously support 
ZBCPs near the transition. These features are proper to 
STIs and distinct from a simple three-dimensional TSC 
mentioned above. Our findin gs s trongly support that the 



observed ZBCPs in Refs.[3Cj, [33j-|35j are originated from 
a helical Majorana fermion in STIs, and they give a firm 
evidence of their topological superconductivity. Our re- 
sults arc summarized in Table |TJ 

First let us briefly review basic properties of parent 
topological insulators. For concreteness, we consider the 
following k ■ p Hamiltonian to describe the topological 
insulators, 

H T i(k) = ma x + v z k z a y + va z (k x s y - k y s x ), 

m = m + mik 2 z + m 2 {k x + ky), (mim 2 > 0), (1) 

where and are the Pauli matrices in spin and 
orbital spaces, respectively. In addition to the time- 
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STI 


BW 


gap 


full 


nodal 


iso 


SABS 


cone/caldera 


ridge/ valley 


cone 


conductance 


DP/ZBP 


ZBP 


DP 



(a) cone/caldera 



TABLE I. Momentum-independent odd-parity paring sym- 
metries in STI. As a comparison, pairing symmetry in BW 
phase of superfluid 3 He is shown. The energy spectrum has 
full gap, nodal or isotropic (iso) full gap. In cases of low and 
intermediate transmissivity of normal metal/STI junctions, 
the line shapes of tunneling conductances show double peak 
(DP) and zero bias peak (ZBP), respectively, (see Fig[3]and 
the corresponding discussions in the text). 



reversal symmetry, we have assumed the mirror symme- 
try MiH T iM\ = ffrik-*-** with Mi = s i: (i = x,y) 
and the inversion symmetry. Although in the above 
is axial symmetric along the z-axis, even if one adds 
higher order terms of fcj (i = x, y) like the warping terms 
471 ] . our results do not change qualitatively. The topo- 
logical phase of this system is classified by the Z2 in- 
variant, (— 1)" = sgn(moTOi), and when Z2 non-trivial 
(momi < 0), the system becomes a topological insula- 
tor. On the surface perpendicular to z-axis, it supports 
the topologically protected Dirac fermion. 

Now consider the corresponding STIs. The STIs are 
described by the Bogoliubov-de Gennes (BdG) Hamilto- 
nian in the Nambu representation (ipaf: VV-I-j — V'lj.i V^t)) 



ffsTi(fe) = (H T i(k) - m)t z + At, 



(2) 



where [i is the chemical potential, t m is the Pauli matrix 
in the particle-hole (Nambu) space, and A is a 4 x 4 ma- 
trix denoting the gap function. For simplicity, we assume 
that A is a constant matrix, which is generally realized 
when the pairing interaction is short-range and attrac- 
tive. Because of the Fermi-Dirac statistics of electrons, 
A satisfies A = s y A s y , thus there are six independent 
pairings, (A, Aa x , Aa z , Aa y s x , Aa y s y , Aa y s z ) (A is in- 
dependent of k). For each independent pairing, we con- 
sider SABS on the surface normal to the z-axis. 

In order to solve the SABS, we consider the semi- 
infinite STI (z > 0) with a flat surface at z = 0. The 
wave function in this system is given by 



V>sti(^ > 0) = }JlUl 



iqjz ik x x ik^y 



(3) 



where qi (I = 1, ••• ,8) is a solution of E = 
E(k x ,k y ,qj) with E(k) being an eigenvalue of Eq.flJJ, 
and ui(k x ,k y ,qi) is the corresponding eigenvector. 
Among the eigenvectors, ipsTi{z) consists of those with 
E(k x ,k y ,qi)/dqi > or Imqi > 0, where the former 
denotes up-going states and the latter denotes localized 
states in the vicinity of z = 0. Postulating the bound- 
ary condition as ipsTi(z = 0) = 0, we can determine the 
coefficients tj and obtain the SABS. 




(d) ridge/valley 
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FIG. 1. (color online) Polar plots of the bulk superconduct- 
ing gap E g for full (a) and nodal (b) gaps, where cone/caldera 
and ridge/valley SABSs are realized, respectively. It is not 
plotted in a certain region for the cases (b), for visibility. 



We find that among the six pairings mentioned above, 
only the three (Aa y s x , Aa y s y , Aa y s z ) support well- 
localized gapless SABSs 011 the surface at z = 0. 
We notice that all of them are odd-parity pairings, 
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-Aa v Sn, 



(P = a x ), and the existence 



following, we focus on Aa y s 
because the result of Aa y s x 
and k y in that of Aa y s y . 



of the SABSs is consistent with the Fermi surface cri- 
terion for odd parity TSCs [2^|. Furthermore, they all 
are odd under (at least) one of the mirror symmetries, 
AiiAayS^Ail = —AaySft for i ^ /a. As illustrated in 
Fig[TJ Aa y s x and Aa y s y have point nodes in the bulk 
spectrum on the the k y - and fc^-axes, respectively, while 
Aa y s z is full gapped. The point nodes change a quali- 
tative structure of the SABSs as is shown below. In the 

(= Af) and A<j y s y {= A n ) 
is obtained by exchanging 

k x 

The obtained SABSs in the STI are illustrated in FigH 
The SABSs appear when m§ < /1 2 . An important fea- 
ture of the SABSs is that there exists a structural transi- 
tion of the energy dispersion. Combined with the nodal 
structure mentioned above, the transition results in a va- 
riety of Majorana fermions, which we call (a) cone, (b) 
caldera, (c) ridge, and (d) valley: For the full gapped 
pairing A = Af, we find that the cone and the caldera 
are possible. For larger values of p, and mi, the energy 
spectrum of the SABS is an axial symmetric and mono- 
tonic function of k[= (fc 2 . + fc 2 ,) 1 / 2 ], and its shape is a de- 
formed cone [Fig[2la)] including higher order terms of k. 
For smaller values of /x and mi, however, a second cross- 
ing of the zero energy appears at finite k and a caldera 
SABS is realized [Fig|2fb)1. This result is consistent with 
that of Refs. 2^, 31, 32| . On the other hand, for the nodal 
pairing A n , we obtain the ridge [Figj2fc)] and the valley 
[FigGJd)], instead. Although the structural transition 
occurs on the same critical line, Majorana fermions in 
this case have the flat dispersion due to the existence of 
bulk point nodes. As a result, the cone (caldera) is de- 
formed into the ridge (valley). We can also show that the 
flat dispersion between the point nodes has a topological 
origin, thus is not accidental [H, \^ . 

Now we show that the structural transition is intrinsic 
to the STI. Due to an argument based on symmetry given 
below, we find that the STI may have a remnant of the 
surface Dirac fermion in the parent topological insulator, 
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FIG. 2. (color online) Structural transition of SABSs in 
STIs. Cone (ridge) and caldera (valley) SABSs are shown in 
the insets. The position of the circle (square) symbol corre- 
sponds to the parameters used in the calculations of tunneling 
conductances with mi = 20.18 eVA (mi = 5.66 eVA), where 
a cone (caldera) or ridge (valley)-SABS is realized. 



and this is the origin of the structural transition. To sec 
this, consider how the superconductivity of A s= f ;I1 may 
affect on the Dirac fermion. When jj, is small and in the 
bottom of the bulk band, the surface Dirac fermion near 
the Fermi energy is well-separated from the bulk band. 
Thus, it can be treated apart, and the problem reduces 
to constructing a pairing term of the Dirac fermions that 
is consistent with symmetry of A g . In particular, the 
induced pairing should have odd mirror parity as A g , 
because the mirror symmetry M.i is a good symmetry on 
the surface at z = 0. However, one find that no pair- 
ing term is allowed to be consistent with the symmetries. 
This means that the Dirac fermion remains to be gapless 
near the Fermi energy when adding A„ , in contrast to the 
case of conventional s-wave pairing [l5|, [5(| . By hybridiz- 
ing with the Majorana cone (ridge) specific to TSCs, the 
gapless Dirac fermion results in a complex caldera (val- 
ley) structure of the Majorana fermions. Now consider 
tuning fi deep into the bulk band. As one increases fi, the 
surface state near the Fermi energy merges into the bulk 
band, and finally disappears. One now obtains conven- 
tional Majorana cone or Majorana ridge because of no 
hybridization of the Dirac fermion. Therefore, a struc- 
tural transition of the Majorana fermions must occur be- 
tween these two limits. 

We note that when the transition occurs, the veloc- 
ity of the Majorana fermions at {k x , k y ) = changes its 
sign. The velocity along the x-direction v is given by 
v = vaA/mo with 



1 - yj\ + 4mi + 4m 2 /} 2 
2mi// 2 



(4) 



line determined by a = is given by jl 2 = l/(— rhi), 
which is shown in Figj2] Only for the case with tooTOi < 
0, the value of a can become zero, namely, topological 
insulator triggers the structural transition of SABS. 

Now we calculate the tunneling conductance of nor- 
mal metal/STI junction, generalizing theories of the tun- 
neling spectroscopy of conventional 5l| and unconven- 
tional [H, HH superconductors. We suppose a free elec- 
tron in the normal metal with the Hamiltonian -ffjv(fc) = 
[(k x + ky + k z )/(2m c ) — hn]ctoSot z . The wave function 
in the normal metal (z < 0) is given by 



ip N (z < 0) 

+E( a ' 



A<TSC L ' 



i-hi 



,(5) 



where Xosr is the eigenvector of H^(k) with orbital 
a and spin s for electron (r = e) or hole (r = h), 
and k cz = \Jk 2 — k 2 = 
khz = y / 2m c (fi N - E) 



k e cos 9, k c = y^meOujv + E), 
k 2 , and a asa ' s i (b asiJ , s ,) is the 
Andrecv (normal) reflection coefficient. The first term 
of the wave function denotes an injected electron, and 
the second (third) one denotes a reflected hole (elec- 
tron) with reflection coefficient a ast7 i s (6 CTSCT ' S ). On the 
other hand, the wave function in the STI side (z > 0) 
is given by Eq. ([3]) with the transmission coefficient 
tj. These wave functions are connected at the inter- 
face (z = 0) by the condition [53j, Vn(0) = V'sti(O) 
and ■unV'nIO) = VstiV'Sti(O), with the velocity operator 



U N(STI) 



dH 



N(STI) 



/dk z 



The above equation 
and tj. Finally, 
the normalized charge conductance G is given by 



determines the coefficients a asr7 i s , b asa > s 



G 



/ d<t> / d8sm26T as (8,(t),eV) 



J0 



E 



2tx r ir/2 



JO 



; (6) 



d0sin20T (TS (0,0,O)|A=o 



where m\ = momi/v 2 and jl = fx/mo- The transition 



with the angle resolved transmissivity T aa (9,<j),E) = 
1 + E ff ' s ' G°Wr's' | 2 ~ \bascr's'\ 2 ) with k x = k COS (j), ky = 
fcsin0, where the energy E of the injected electron is 
fixed at the bias voltage eV. 

In the following, the band mass of the normal metal 
is fixed as m TO 2 = 1 f° r simplicity, and we set A 
().() meV and mi = —0.59 or rh\ = —0.17. The other 
parameters are the same as those used in Ref. [3(|, i.e., 
too = -0.28 cV, m 2 = 56.6 eVA 2 , v z = 3.09 eVA, v = 4.1 
eVA and ft = —1.8. We control the transmissivity of the 
normal metal/STI interface by changing the value of /in- 
The transmissivity becomes maximum in this model for 
Mjv/a* ~ 0-6 since the magnitude of Fermi momentum 
in the normal metal coincides with that in STI. As /xjv 
increases, the magnitude of transmissivity is reduced. 

The obtained tunneling conductances G/Gn near the 
structural transition of SABSs as functions of bias volt- 
age eV/A are shown in FigEI With the decrease of the 
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FIG. 3. (color online) The normalized tunneling conductances 
G/Gn near the structural transition of SABSs as functions 
of bias voltage eV/A for the cone, caldera, ridge and valley 
SABSs. The values denoted in the panel are of \i/ \xn for each 
line. 



magnitude of transmissivity (/ijv//-* = 60, 1200), the ro- 
bust ZBCPs appear stemming from the gapless SABSs 
in Figl^l Only for the low transmissivity case with 
Hn/h = 1200 as shown in Figj3](a), G/Gn has a double- 
peak structure. The latter is consistent with the fact that 
the corresponding surface local density of states does not 
have a zero energy peak but double peak structure [3fj| . 
In junctions with high transmissivity with hn / H = 0.6, 
we obtain G/Gn ~ 2 for \eV\ ~ which is also consistent 
with the fact that an injected electron is almost perfectly 
reflected as a hole due to Andrccv reflection. 

We now focus on STI with Af. It is noted that the 
difference of the line shapes of G/Gn between Fig(3Ja) 
and Fig|3jb) can be understood from the different types 
of SABSs. In the case of FigJ3Jb) , a caldera-SABS is re- 
alized, as shown in FigJ^b). From Eq.Q, the slope of 
the dispersion of the SABSs at k — becomes gradual 
near the structural transition. This enhances the surface 
local density of states at E = and makes the ZBCP 
for the tunneling conductance in the STIs. Thus the 
G/Gn at eV = is enhanced in comparison with that 
for the cone-shaped SABS. As a result, even in the low 
transparent limit /xjv/m = 1200, no double-peak struc- 
ture of G/Gn appears in Figl^b). The present feature 
is different from preexisting 3D TSCs with spin-triplet 
p-wave pairing realized in BW phase of supcrfluid 3 He, 
where in contrast to Figs. 2(a) and (b), the energy dis- 
persion of the SABS becomes a conventional Majorana 
cone [4ll - [46| . In this case, the angle resolved transmis- 




FIG. 4. (color online) The normalized tunneling conductance 
G/Gn as a function of bias voltage for BW state. 



sivity T{6, 0, eV) is given by 



T(0,<j>,eV) 



~~2 



E 

a=±l 



l + a N \T\ 2 + {aN-l)\T\ A 
|l + (l-cr A r)r 2 exp(-2i6»s) 



(7) 

with the transmissivity at the interface <jn given by 
<7 N = 4 cos 2 (9 / (4 cos 2 6> + Z 2 ) [12] and T = A/(eV + 
y/ (eV) 2 — A 2 ). Z is a dimensionless constant that con- 
trols cat, and A is the superconducting gap in this sys- 
tem. The resulting tunneling conductance never shows a 
ZBCP [Io| as shown in Figf4j This difference comes from 
the absence of the structural transition. 

Next, we consider STI with A n , where the resulting 
SABS has a quasi-one dimensional energy dispersion. In 
the x-direction, SABS has a flat dispersion as mentioned 
before [Figl^c)]. The present flat dispersion of the SABS 
makes a ZBCP in G/Gn for arbitrary lower transmissiv- 
ity, as shown in Figj3jc). When a valley-cone is realized 
as the SABS [FigHfd)]. G/G N at eV = is enhanced 
[FigHd)]. 

Finally, we compare our results with the experimen- 
tally observed tunneling spectroscopy in Cu x Bi 2 Sc3. The 
tunneling conductance in Au/Ag/Cuo.3Bi2Se3 junction 
has been observed in Refj30| . From the lattice constants 
of Au and Ag (a ~ 4A) [54[ , the Fermi momentum of the 
normal metal is estimated as kp ~ n/a ~ 1A _1 , which 
corresponds to [In / ^ ~ 100 in our model. While in the 
actual system, a barrier layer suppressing transmissivity 
could be formed between normal metal and STI, it can 
be taken into account as an effective increase of /zjv//i. 
Therefore, the experimental result in Ref. [13] should be 
compared with ours with > 100. From Fig|3l we 

find that the experimentally observed ZBCP is consistent 
with Af and A n , both of which support ZBCPs originat- 
ing from Majorana fermions on the normal metal/STI 
interface. 

In conclusion, we have developed a theory of the tun- 
neling spectroscopy of STL We have clarified the struc- 
tural transition of the energy dispersion of the SABS, i.e., 
cone-caldera and ridge-valley transitions, which stems 
from remaining metallic surface states of parent topo- 
logical insulator. In the vicinity of the structural tran- 
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sition of SABSs, even in the full-gap superconducting 
case, the line shapes of tunneling conductance show ro- 
bust ZBCPs. On the other hand, a typical 3D topological 
superconductor with pair potential realized in BW phase 
in superfluid 3 He, never shows a ZBCP. Our obtained 
results serve as a guide to explore novel topological su- 
perconductors with Major ana fermions [HH - 57 1 . 
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